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ABSTRACT 

Color-magnitude diagrams (CMDs) derived from Hubble Space Telescope 
(HST) Advanced Camera for Surveys F606VI^, F814VI^ photometry of 6 globular 
clusters (CCs) are presented. The six CCs form two loose groupings in Calac- 
tocentric distance (Rqc): IC 4499, NCC 6426, and Ruprecht 106 at -15-20 kpc 
and NCC 7006, Palomar 15, and Pyxis at — 40 kpc. The CMDs allow the ages to 
be estimated from the main sequence turnoff in every case. In addition, the age 
of Palomar 5 (Rqc ^ 18 kpc) is estimated using archival HST Wide Field Plane- 
tary Camera 2V^I photometry. The age analysis reveals the following: IC 4499, 
Ruprecht 106, and Pyxis are 1-2 Cyr younger than inner halo CCs with similar 
metallicities; NCC 7006 and Palomar 5 are marginally younger than their inner 
halo counterparts; NCC 6426 and Palomar 15, the two most metal-poor CCs in 
the sample, are coeval with all the other metal-poor CCs within the uncertainties. 
Combined with our previous efforts, the current sample provides strong evidence 
that the Calactic CC age-metallicity relation consists of two distinct branches. 
One suggests a rapid chemical enrichment in the inner Calaxy while the other 
suggests prolonged CC formation in the outer halo. The latter is consistent with 
the outer halo CCs forming in dwarf galaxies and later being accreted by the 
Milky Way. 
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Subject headings: globular clusters: general — Galaxy: formation 



Introduction 



The Galactic globular clusters (GCs) have long been recognized as useful probes of the 
formation and evolution of the Galaxy. Their spatial coherence allows their locations in age, 
metallicity, and Galactocentric distance (Rqc) space to be established with some accuracy. 
Their numbers allow them to be used as probes of the formation timescale and chemical 
enrichment history of the Milky Way halo, thick disk, and central bulge. 

One of the earli est papers to examine the age-metallicity relation (AMR) of the Milky 
Way GC system was ISarajedini fc King ( 119891 ). They calculated ages for 32 GCs using the 
magnitude difference betwe en the horizont al branch (HB) and the m ain sequence turnoff 
(MSTO). Taking a cue from iGrattonI fll985l ). ISarajedini fc Kin3 fll989h compared the AMR 
for GCs inside 15 kpc of the Galactic center with that of GCs outside 15 kpc. They found a 
statistically significant difference between the two: the AMR of GCs with Galactocentric dis- 
tance (Rqc) < 15 kpc was shallower t han that of GC s wit h Ran > 15 kpc. ISaraiedini fc King 
( 1989 ) argued, based on the work of Larson ( 1972 ) and Tinsley fc Larson ( 19781 ). that the 
inner halo GCs formed in a region of higher gas and dust density, thereby accelerating GC 
formation and chemical enrichment. In contrast, the outer halo GCs formed in lower density 
sub-galactic fragments in a slower, more chaotic f ashion akin to the fragmentation and accre- 
tion scenario advocated by ISearle fc ZinnI ( 119781 see also Navarro et al. 1997) and featured 
in models of g alaxy formation based on the cold dark matter paradigm with a cosmoloRical 
constant (e.g., Izentner fc Bullock 2003 : Robertson et al. 2005 : Font et al. 2006 ). 



Subsequent work bv lChabover et al I (|l992r ). ISaraiedini. Lee, fc Led (|l995r> . and lChaboyer et al. 



I 1 j , ^ ^^^^^^^^^^^^^^ ^ . ^ , — ^^^^^-.-.-^^-^^^^^^g^^^— — — ,^ ^ , ^, 

( 119961 ) reinforced the results of lCrattonl ( 119851 ) and ISarajedini fc KingI ( 119891 ) using larger and 



more reliable data sets of GC ages as well as a variety of techniques for measuring them. The 



most recen t studies to deterni ine ages from a large sample of GCs, by iMarm-Franch et al. 



(l2009h and botter eTaP (l201oh . rely upon the uniform photometric data set provided by the 
Hubble Space Telescope (HST ) Advanced Camera for Surveys (ACS) Galactic GC Treasury 
project ( ISarajedini et al.ll2007l ). T his data set provides d eep photometry for 65 GCs of up to 
^7 magnitudes below the MSTO ( lAnderson et al.ll2008l ). Based on a careful analysis of the 



Based on observations with the NASA/ESA Hubble Space Telescope^ obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555, under program 
GO-11586. 
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ages and metallicities of the Galactic GC Treasury project GCs, iMarm-Franch et al.l (120091 ) 
showed the clearest evidence thus far that the G C AMR splits into two distinct branches. 
This result was reinforced by iDotter et al.l ( 120101 ) , who added the six most distant (known) 
Galactic GCs to the GC Treasury data set and used the combined sample to study the cor- 
rel ations between HB morphology and a variety of GC parameters. The analysis presented 
by iForbes fc Bridges! (120101 ) using data assembled from the literature for 93 GCs largely 
reinforces these conclusions. 

The HST/ACS Galactic GC Treasury project preferentially targeted nearby GCs in or- 
der to maximize the depth that could be achieved with one orbit per filter. One consequence 
of this strategy is that it left out much of the outer halo. The most distant GC observed as 
part of the Treasury project is NGC 4147 at Rqc ^ 21 kpc. Yet there are many GCs beyond 
this distance that can provide more leverage on the formation and subsequent evolution of 
the outer Galactic halo. 

A number of investigators have used HST Wide Field Planetary Camera 2 (WFPC2) 
photometry to derive the ages of GCs with Ra n > 50 kpc. [Harris et al.l ( 119971 ) presented 
deep photometry of NGC 2419 (Rqc ^ 90 kpc: lHarrisl[l996l . 2010 revision) and concluded 
that it is coeval with the similarly metal- poor GC M92 based on a differential comparison of 
their color-magnitude diagrams (CMDs). IStetson et al.l ( 119991 ) obtained deep photometry of 
Palomar 3, Palomar 4, and Eridanus (all with Rqc > 90 kpc) and concluded that the three 
are ^1.5-2 Gyr younger than the relatively nearby GCs M3 and M5. Their conclusion rests 
on the assumption that the outer halo GCs have similar chemical abundances ([Fe/H] and 



a/Fe]) to M 3 and M 5 ; this assumpt i on ha s since been bo r ne out for Palomar 3 (i Koch et al. 



20091) and Palomar 4 ( iKoch fc Cotd l2010l ). iDotter et al.l (l2008l ) presented CMDs of Palo- 
mar 14 (Rqc ^ 75 kpc) and AM-1 (Rqc ^ 125 k pc) and concluded th at both are 1.5-2 Gyr 
younger than M3, in keeping with the results of IStetson et al.l ( 119991 ). For the purposes of 
the present study, the 6 known GCs with Rqc > 50 kpc are sufficiently well studied to firmly 
mark their place in the GC AMR. Chemical abundance information remains rather sparse 
for the most distant GCs, but what is known is consistent with what has been inferred from 
the CMDs. 

Thus far then, deep and fairly homogeneous HST CMDs in the V and / (or equivalent) 
filters exist for the majority of GCs with Rqc ^ 15 kpc and all those (known) beyond 50 kpc, 
allowing precise and internally consistent ages to be measured. What remains is to target 
the GCs within 15 < Rqc ^ 50 kpc so that the complete radial extent of the Galactic GCs, 
and what it can tell us about the formation of the Galaxy, can be probed. The present study 
is the first major step toward fulfilling this goal, wherein we present deep, homogeneous HST 
ACS photometry, and use the resulting CMDs to estimate the ages, of six outer halo GCs: 
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Pyxis, Ruprecht 106, IC 4499, NGC 6 426, NGC 7006, an d Palomar 15. Of the six, only 
IC 4499 has previous HST photometry (jPiotto et al.ll2002l ). and it is not of sufficient depth 
to aUow an accurate age determination. In addition, a new a ge determination for Palom ar 5 
(Rgc ^ 18 kpc) using archival HST/WFPC2 VJ data from ICrillmair fc SmithI (120011) has 
been performed and included in the subsequent discussion. 

The remainder of this paper proceeds as follows: the observations and data reduction 
are described in §2, the CMDs are presented in §3, the available information on the chemical 
abundances, distances, and reddening values of the sample is reviewed in §4, new age deter- 
minations based on isochrone fitting are presented in §5, the Galactic GC AMR is updated 
and compared with theoretical work in the context of the formation of the Galaxy in §6, 
§7 provides some discussion of the results, and §8 summarizes the findings presented in the 
paper. 



Observations and Data Reduction 



The observations were obtained with the HST/ ACS Wide Field Channel (WFC) under 
program number GO-11586 (PI: Dotter) during Cycle 17. Table [1] shows the observing log. 
The ACS images were retrieved from the HST archive and calibrated using the pipeline bias 
and fiat-field procedures. All clusters, except for Palomar 15, were imaged over two orbits 
with the first orbit devoted to F606VI^ and the second to F814:W. Palomar 15 is the most 
distant and reddened in the sample and it received a total of five orbits: two in F606VI^ and 
three in F8UW. 

Photometry and astrometry were extracted from the images using the programs de- 
veloped for GO-1077 5, the ACS Survey of Galactic Globular Clusters, as described by 
Anderson et al.l ( 120081 ). The observing strategy in the current program was designed to fol- 



Table 1. GO-11586 Observations 



Cluster 



Data set 



Date 



RA 



Dec 



PA_V3 



F606W 



F814W 



Pyxis 

Ruprecht 106 
IC 4499 
NGC 6426 
NGC 7006 
Palomar 15 



jbl601 
jbl602 
jbl603 
jbl604 
jbl605 
jbl606 



11/10/2009 
04/07/2010 
01/07/2010 
04/08/2009 
05/10/2009 
16/10/2009 



09h07m57.7s 
12h38m40.1s 
151i00ml8.5s 
17h44m54.6s 
21h01m29.4s 
161i59m50.9s 



-37:13:17.0 

-51:09:00.9 

-82:12:49.5 

+03:10:13.0 

+16:11:14.4 

-00:32:17.9 



76.86 

287.34 

235.97 

294.67 

275.03 

261.25 



50s, 4 X 517s 
55s, 4 X 550s 
60s, 4 X 603s 
45s, 4 X 500s 
45s, 4 X 505s 
10s, 4 X 550s 
65s, 4 X 550s 



55s, 
60s, 
65s, 
50s, 
50s, 
10s, 
25s, 
55s, 



X 557s 
X 585s 
X 636s 
X 540s 
X 545s 
X 500s 
X 560s 
X 525s 
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l ow th e identical pattern of exposure times and dithers as GO-10775. Since lAnderson et al. 
( 120081 ) give an extremely detailed description of the data reduction process, only those as- 
pects that have changed in the interim will be mentioned here. 

Before performing photometry on the individua l images, they were correc ted for charge 
transfer efficiency using the algorithm developed by lAnderson fc BedinI ( 120101 ) . The instru 



menta l photometry was then calibrated to the HST VEGAmag sys tem following iBedin et al. 



( 2005 ) and using the aperture corre ctions from lSirianni et al.l ( 120051 ) as described bv ^Sarajedini et al 
(120071) and lAnderson et al.l (120081 ). The VEGAmag photometric zeropoints were obtained 
from iBohlinI ( 120071 ). The photometric catalogs and supporting data files from the reduction 
process of the 6 GCs presented herein will be made publicly available through the same 
archive that will host the ACS GC Treasury database. 



3. Color-Magnitude Diagrams 

The F606W, F606W - F8UW CMDs are presented in Figures □ through El Typical 
photometric errors are demonstrated towards the left side of each figure. In cases where 
ground-based I data are available, we also make a direct comparis on with the ACS dat a 



converted to V and / using the empirical transformations provided by lSirianni et al.l ( 120051 ). 
The ground-based comparisons indicate, in all cases, that the HST photometric system 
maintains a high degree of homogeneity post-SM4 and that the transformations to ground- 
based V and / magnitudes perform well. 

The CMDs reveal several striking features. In particular, it can be seen that Ruprecht 106 
and IC 4499 have strong binary and blue straggler sequences. NGC 6426, Palomar 15, and 
Pyxis show signs of differential reddening and, indeed, these are the most reddened clusters 
in the present sample. Palomar 15 and Pyxis are sparsely populated clusters; the red giant 
branch (RGB) of Pyxis loses coherence for F606VI^ < 19. 



4. Previous studies 



This section provides a brief review of noteworthy studies targeting the GCs considered 
in the present study. A few basic parameters are collected in Table El 
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Fig. 1.— Left: The F606Vr-F814Vr CMD of Ruprecht 106. Right: T he ACS data f crosses) 
conver ted to V and /, compared with the ground-based V and / data of lSarajedini fc Lay den 
(119971 . open circles). 



-7- 




Fig. 2.— Left: The F606W - F8UW CMD of IC 4 499. Right: The AC S data (crosses), 
converted to V and /, compared with the V, I data of I Walker et al.l (120111 . open circles). 
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Fig. 4.— Left: The F606W - F8UW CMD of NGC 642 6. Right: The ACS data (cr osses), 
converted to V and /, compared with the V, I data of iHatzidimitriou et alJ (119991 . open 
circles) . 
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F606W-F814W V-I 



Fig. 5.— Left: The Fm6W - F8UW CMD of NGC 7006. Ri g ht: T he ACS data (crosses), 
converted to V and /, compared with the V, I data of IStetsonI (l200fll . open circles). 
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Fig. 6.— The FQOQW - F8UW CMD of Palomar 15. 



Table 2. Literature values of basic cluster parameters 



Cluster 


DMv 
(Harris) 


E{B - V) E{B - V) 
(Harris) (SFD98) 


[Fe/H] 
(Harris) 


[a /Ye] 


IC 4499 


17.09 


0.23 


0.224 


-1.6 




Ruprecht 106 


17.25 


0.20 


0.174 


-1.5'' 


- 0.0" 


NGC 6426 


17.70 


0.36 


0.346 


-2.26 




NGC 7006 


18.24 


0.05 


0.082 


-1.55^ 


+0.23^ 


Pyxis 


18.65 


0.21 


0.327 


-1.4^ 




Palomar 15 


19.49 


0.40 


0.394 


-2.0^^ 




Palomar 5 


16.92 


0.03 


0.056 


-1.3^= 


+0.16" 



References. — Sources of data a re listed in th e column headers except 
where otherwise specified. Harris- lHarris fll996h : S FD98-ISchlegel et al. 



(Il998h: ral P. Fr ancois (2010) priv. comm . : ( iKraft et al I (ll 9m : 



Smith et all (12002) 



ic^ IPalma et all J200oh : {d) IPa Costa fc Armandrofj (Il995h : (e) 
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4.1. NGC 6426 



The first modern CMD of this cluster was presented by lZinn fc BarnesI (119961 ) based on 
imaging with the 1-m telescope at Cerro Tololo Inter- American Observatory. Their CMD 
exten ded ^3 mags below t he HB and confirmed the id e ntity of 12 variable stars in the clus- 



ter. 



Hatzidimitriou et al.l ( 119991 ) and IPapadakis et al.l ( 12000 1) published BVRI photometry 



extending to about 1 mag below the MSTO and an ex tensive study of its RR Lyrae variable 
population, respectively. iHatzidimitriou et al.l ( 119991 ) concluded that NGC 6426 is coeval 
with the mean age of the metal-poor Galactic GCs. 



4.2. IC 4499 



The first CCD photometry of IC 4499 was published by ISarajedinil (119931 ) wherein a 
B — V CMD extending ~3 mags bel ow the HB wa s prese nted. IC 4499 is exceptional for 
its prodigious RR Lyrae population. ISuntzeff et al.l ( ll99ll ) compiled luminosity-normalized 
RR Lyrae populations for 77 Galactic GCsEI On a scale where N^^^^ = 0.3 for 47 Tuc and 
N^^ = 56.4 for M3, the N^^^^ value of IC 4499 is 88.7. The only GC with a higher value of 

is Palomar 13, but with only four RR Lyrae stars its specific frequency is likely to be 
infiuenced by small number statistics. 

A comprehensiv e study of the CMD and RR Lyrae population of IC 4499 was performed 
by IWalker fc Nemed (119961 ). who derived a mean reddening of E(B— V)=0.22 ± 0.02 based 
on four independent estimates and [Fe/H]= —1.65 ± 0.10 on the IZinn fc WestI (Il984r) scale. 
Adopting a reddening and metallicity close to these values, iFerraro et al.l ( 119951 ) conclude 
that IC 4499 is younger (by 3-4 Gyr) than most globular clusters at its metallicity. 



Kunder et al.l ( 120111 ) presented an updated study of the RR Lyraes in IC 4499 and found 



that thei r period chan g e rate s were an order of magnitude larger than predicted by stellar 
models. I Walker et al.l ( 120111 ) presented multi-band, ground -based photometry o f IC 4499. 
Their CMD reaches ^ 3 mag below the MSTO from which I Walker et al.l ( 120111 ) estimated 
an age of 12 ± 1 Gyr. 



The quantity used by ISuntzeff et al. I (|l99l[ ). ^RR^ is an integrated luminosity- averaged quantity: the 
number of RR Lyraes a GC would have if the GC itself had My = —7.5. 
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4.3. Ruprecht 106 



Discovered b y Ruprecht on photographic plates as he was searching for open clusters 



dA 



ter et al 



19611). Rupre cht 106 existed in relative obscurity until the first CCD photometry 



by lBuonanno et al.l (Il990l ) showed that it is metal-poor with a red HB and, t hus, was likely to 
be you nger than other GCs with similar metallicities. Subsequent work by iBuonanno et al. 
( 1 19931 ) corroborated this result. 



The metal abunda nce of Ruprecht 106 has been a source of controversy since the work of 
Buonanno et al.l ( jl990l ) . The abundance studies performed thus far have fallen into two broad 
categories, those that find a metal abundance of [Fe/H] ^ —1.9 ( IBuonanno et al.lll990l . Il993l : 



Sarajedini fc LaydenI 119971 ) base d on photometric indicators and those that c onclude [Fe/H] 



— 1.6 based on spectroscopy ( IDa Costa et al. 



1992 



Francois et al.l 119971 ). iKaluzny et al. 



( 119951 ) found that Ruprecht 106 and M3 have similar metallicities based on time series 
photometry of 12 RR Lyrae variables and, separately, that [Fe/H] > —1.6 based on the 
relative locations of the RGB bump and the HB. 

The discrepancy in metallicity estimates may be due to t he unusually low [a/Fe] rati o 
of Ruprecht 1 06 co mpared to other metal-poor Galactic GCs ( jPritzl. Venn, fc IrwinI l2005l ) . 
Brown et al.l ( 119971 ) claimed [Fe/H] =2::^ —1.45 and [0/Fe] ^ based on high-resolution spec- 
tra of 2 red giants. VLT/UVES spectroscopy of 6 red giants indicates —1.5 < [Fe/H] < —1.45 
and —0.1 < [a/Fe] < +0.1 (P. Francois, 2010, private communication). 



4.4. Palomar 15 



The e arliest CCD p hotometry of Palomar 15 dates back to the work of ISeitzer fc Carney 



( 11 9901 ) and iHarrisI ( jl99ll ). CMDs presented in both of these studies extend from the tip of the 
RGB to about 3 mags below the HB and reveal an HB morphology that is blue- ward of the 
RR Lyrae instability strip. In addition, both studies conclude that Palo mar 15 suffers from 



ajiigher-than-expected amount of line-of-sight reddening. Whereas the 
(Il982h reddening maps predict E(B— V) ^0.1, the CMDs of ISeitzer fc Carneyl 
HarrisI ( 1991r) suggest a redd ening closer to E(B— V) ^ 0.4, consistent with the reddening 



Burstein fc Heiles 
( Il990h and 



maps of lSchlegel et al.l ( 119981 ). Given this value, photometric indicat ors predict [Fe/H] ^ -1.9. 
The o nly spectroscopic value of the metallicity was published by iDa Costa fc Armandrofi 
( 119951 ). who found [Fe/H]= —2.00 ± 0.08 based on the strength of its Calcium triplet lines. 
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4.5. NGC 7006 



NGC 7006 was recog nized early-on as a prime example of the second parameter effect. 
Sandage fc WildevI fll967l ) noted that that M13, M3, and NGC 7006 aU have very similar 
metallicities ([Fe/H] ^ —1.6) but exhibit a range of HB morphologies, with M 13 displaying 
the bluest HB and NGC 7006 th e reddest. Despit e its im portance in this regard, the deepest 
published C MP remains that o f iBuonanno et al.l ( ll99ll ). which extends only ^1 mag below 

1998|) estimated [Fe/H]= —1.55 based on high-resolution spectra 



the MSTO. 



Kraft et al. 



of 6 stars. iKirby et al.l ( 120081 ) found similar a similar average value from medium-resolution 
spectra of 20 stars: [Fe/H]= -1.59 ± 0.03. 



4.6. Pyxis 



Discovered by lWeinbergerl ( 119951 ) while searching; t he optical sky surveys, the Pyxis glob- 



ular c luster was confirmed as a GC by Da Costal ( 119951 ) and, soon thereafter, by llrwin et al. 



( 119951 ) and ISarajedini fc Geislerl ( 119961 ) . These studies rev ealed a sparsely populated CMD 
indicative of an intermediate metallicity GC with a red HB. iSaraiedini fc Geisler used 
the simultaneous reddening and metallicity method of ISarajedini fll994h to derive [Fe/H] = 
— 1 .20 =b 0.15 an d E(B — V) = 0.21 ± 0.03 for Pyxis. Spectroscopy of 6 bright giants in Pyxis 
by IPalma et al.l ( 120001 ) yields a mean metallicity of [Fe/H] = — 1.4 ± 0.1, broadly consistent 
with the abundance derived by ISarajedini &: Geislerl ( 119961 ) . 



4.7. Palomar 5 



The first CCD-based photometry of Palomar 5 was published by ISmith et al.l ( 119861 ) 
in the B and V filters. As is typical for the Palomar clusters, the CMD of Palomar 5 is 
relatively sparse on the RGB and HB. Des pite the fact that their CMD clearly delineated 
the MSTO of Palomar 5, ISmith et al.l (|l986r) were unable t o deter n iine a reliable age es t imate 
for the cluster using the isochrones of IVandenBerg fc Belli ( 119851 ). ICrillmair fc SmithI ( 120011 ) 
presented a deeper, more precise CMD of Palomar 5 based on HST/WFPC2 observations in 
the F555VI^ and F814VI^ filters, which they converted to ground-based V, I. The resulting 
CMD extends over 7 magnitudes below the MSTO and was used by iGrillmair fc Smith 
(|2001r) to study the mai n sequence luminosity function of Palomar 5. Since the study of 
Grillmair fc SmithI (l2001h. Palomar 5 has been the subject of extensive investigations focused 



on its tidal tails ( IGrillmair fc Dionatos 



2006 



Odenkirchen et al.ll2009r>. Spectroscopy of four 



giant stars in Palomar 5 yields a mean metallicity of [Fe/H] ^ —1.3 ( ISmith et al.ll2002l ) which 
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is consistent with the abundance inferred from the CMD (ISmith et al.lll986r>. The reddening 
of Palomar 5 is thought to be relatively low, as noted by Smith et al.l ( 1986 ). who quote a 
value of E{B -V) = 0.03. The maps of ISchlegel eTaP (ll998|) give a reddening of E{B - V) 
= 0.056 in the direction of Palomar 5. 



5. Age determinations from isochrone fitting 



Potter et al.l ( 120101 ) presented isochrone fits to the photometric catalog of the ACS Sur- 



vey of Galactic GCs, excluding only those known at the time to possess mult i ple, d istinct 



stellar populations. The fits were performed with isochrones from iDotter et al. 



procedure employed in the isochrone fits presented in this section follows that of 



(l2008h. The 



Potter et al. 



( 120 lOL Sec. 4.2). The distance, reddening, and chemical composition of the isochrones are 
initially set to the best available literature values (see Table [2j) and these are adjusted only 
when necessary to achieve the best simultaneous agreement between the mo dels and the 



CMD on the unevolved main sequence and the base of the red giant branch. I Potter et al. 



( 120 lOl . see their Figures 6 and 7) showed that the distances derived in this manner were 
consistent with the absolute magnitudes of HB stars derived from RR Lyraes with the la 
uncertainties and that the derived metallicities were typically within 0.1-0.2 dex of spectro- 
scopically determined values, particularly in the range —2 < [Fe/H] < —1.5. The values 
derived from isochrone fitting for the present sample are compared to the values obtained 
from the literature at the end of this section. 

Once the isochrones were registered in this manner, the age is estimated by selecting the 
isochrone which most closely matches the shape and position of the subgiant branch. The 
isochrone fits are performed 'by eye' and the quoted uncertainties refiect the width of the 
data in the CMD and any mismatch between the data and the slope of the models through 
this region. As such, these uncertainties are suggestive of goodness-of-fit but lack statistical 



significance. IChaboyer fc KraussI ( 12002 ') report that uncertainties in stellar models lead to 



an error of ±3% in derived ages if a star's properties are known exactly. At an age of 13 Gyr, 
3% corresponds to ^ 0.4 Gyr. Adding this value in quadrature with a typical uncertainty of 
1 Gyr due to the fi t (see Table [3j) and ^ 0.5 Gyr for 0.2 dex uncertainty in metallicity (e.g.. 



Dotter et al.ll2008l ) leads to a realistic error budget of ^ 1.2 Gyr. 



As can be seen in the figures that follow, age effects are only relevant from ^ 2 mag- 
nitudes below to ^ 1 magnitude above the MSTO: greater depth in the CMD provides a 
tighter constraint on the models. The CMDs of Ruprecht 106, IC 4499, and NGC 7006 
provide excellent constraints; differential reddening is a limiting factor in NGC 6426, Palo- 
mar 15, and Pyxis. The sparseness of the Pyxis and Palomar 15 CMDs also reduces the 
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precision of the isochrone fitting. 



Stars brighter than the subgiant branch are saturated in the WFPC2 observations of 
Palomar 5. To enable a more robust age determi nation, t he W FPC2 data have been sup- 
plemented with ground-based observations from Stetson ( 20001 ). For the ^20 stars that 
appear in both catalogs, the average offsets are less than 0.01 mag in V and /; the WFPC2 
data we r e adju sted to match the photometric zeropoints of the standard stars measured by 
StetsonI fcoooh . 



The isochrone fits are shown in Figures [TJ [HI and [9l the results of isochrone fitting are 
presented in Table [31 For comparison purposes, the reddening and distance values listed 
in Table [3] have bee n transformed t o stan dard DMy and E{B — V) using the extinction 
coefl&cients given by ISirianni et al.l ( 120051 ) . The CMDs in Figures [7] and [8] were cleaned 
to improve the clarity of the fits by selecti ng stars accord i ng; to photometric errors and 
other diagnostic information as described by [Anderson et al.l ( l2008[ . Section 7) . As check on 
the results of isochrone fitting, synthetic zero age HB model sequences for the color ranges 
appropriate for each GC are plotted in Figures [7H9] as dashed lines. Generally speaking, these 
comparisons lend confidence to the results of the isochrone fitting. However, note that the 
models appear to inadequately represent those GCs that harbor both red and blue HB stars 
and, further, that the scarcity of HB stars in Pyxis and Palomar 5 make these comparisons 
of relatively little use. 

Comparison of the [Fe/H] values listed in Tables [2] and [3] indicate that the literature 
values and those adopted in the isochrone fits differ, at most, by 0.1 dex in [Fe/H]. The largest 
departures are in cases where no information regarding [a/Fe] is available. Comparison of 
the derived distance moduli shows that 4 of 7 GCs considered here differ by less than 0.05 
mag. The largest discrepancy in distance modulus , that o f Rupr echt 106, also corresponds 
to the largest discrepancy in [Fe/H] between the [Harrid ( jl996[ ) catalog and the adopted 
spectroscopic value. The Harris catalog distances are based a the relation between the 
absolute magnitude of the RR Lyrae stars and the mean [Fe/H] of the GC. 



6. The age-metallicity relation of Galactic globular clusters 



As described in §2, all aspects of the observations and data reduction of the 6 GCs 
observed with ACS in program GO- 11 586 were designed to be homogeneous with the ACS 
Survey of Galactic GCs. The observati ons and data reduc tions of Palomar 5 are c onsis- 
tent with the outer halo GCs studied bv lStetson et al.[ (119991 ) and [Dotter et al.[ ( 120081 ). The 
isochrone analysis of all 7 GCs presented herein is consistent with the approach used in 
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Dotter et all ( 120081 ) and IPotter et all (l2010r>. Hence, the results in Table [3] can be incorpo- 
rated into the AMR assembled by iDotter et al.l (120101 ) . 



The upper-left panel of Figure [TO] shows the AMR of IPotter et al.l (120101 ) with filled 
circles indicating Rqc ^ 8 kpc and diamonds indicating Rqc > 8 kpc. The current sample 
are plotted as open diamonds to highlight where they fall on the AMR. In order to give a 
sense of the uncertainties in the diagram, a 'typical' error bar representing ±0.2 dex in [Fe/H] 
and ±1 Gyr in age is plotted in every panel of Figure [101 The remaining panels compare 
the GC AMR to various model predictions: the upper-right panel shows the si mple closed 



box ch emical evolution model for the Sagittarius dwarf spheroidal galaxy from ISiegel et al 



the lower-left pane l show s the chemical evolution models of the Magellanic Clouds 
from IPagel fc Tautvaisiend (119981): and the lower-r ight panel shows the mean trend from the 
GC formation model of iMuratov fc GnedinI (120101 ). These comparisons are discussed in the 
following paragraphs. 

The closed-box model portrayed in the upper-right panel of Figure [TOl was taken from the 
study of the star formation hist ory of M 54 and the central field of the Sagittarius (hereafter 
Sgr) dwarf spheroidal galaxy by ISiegel et al.l (120071 ) . The star formation history was derived 
with the same stellar models used in this paper and photometry from the ACS Survey of 
Galactic GCs. It has already been demonstrated by iLayden fc Saraiedinil (l2000l) that the 
GCs associated with Sgr follow the same trend as its field stars. iForbes fc Bridged ( 120101 ) 
reiterate that the AMR of Sgr is essentially identical to that of the putative Canis Major 
(CMa) dwarf galaxy and argue that ^ 20 GCs originated with one or the other. The upper- 
right panel of Figure [10] indicates that the outer halo GCs follow an AMR consistent with 
dwarf galaxy chemical evolution. However, it appears that the outer halo GCs formed in a 
variety of environments within which chemical enrichment proceeded at different rates; one 
characteristic of Sgr or CMa is insufficient to explain the observed AMR. 

The GC AMR is compared to the A MRs of the SmaU and Large M ageUanic Cloud (SMC 
and LMC) chemical evolution models bv lPagel fc Tautvaisiend ( 119981 ) in the lower-left panel 
of Fig ure [Tot The LMC AMR is not markedly different from that of M54/Sgr, as pointed 
out by iForbes fc Bridged (120101 ). The SMC AMR tracks the majority of GCs with ages less 
than 12 Gyr better than either the LMC or M54/Sgr. 

The comparisons with Sgr and the MCs indicate that a collection of dwarf galaxies with 
a range of AMRs ar e likely to have co n tribut ed to the Galactic GC popu lation, consistent 
with the findings of IForbes fc Bridged (120101 ). iMuratov fc GnedinI (120101 ) devised a semi- 
analytical model of the formation of a massive host galaxy's GC population and trained their 
model on the Galactic GC metallicity distribution. The model is built on top of cold dark 
matter galaxy formation simulations and assigns a globular cluster formation probability and 
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chemical enrichment history to each sub-halo that is accreted by the host galaxy. Among 
their stated goals is the development of a bimodal color distribution in the GC population 
without requiring two distinct formation mechanisms. While the model achieves this goal 
by design, it also fits the age-metallicity distribution of the outer halo GCs remarkably well. 
The lower-right panel of Figure [10] compares the GC AMR with the mean trend (squares) 
and standard deviation (indicated by e rror bars) of [Fe/H ] for .5 Gyr age bins of nearly 
10,000 simulated GCs drawn from the llVIuratov fc Gnedin ( 2010 ) model. While the outer 
halo GCs are well represented by the mean trend, the model is discrepant with the metal-rich 
GCs in the inner Galaxy, i.e., the bulge and thick disk GCs. Muratov & Gnedin attribute 
this discrepancy to the lack of a metallicity gradient within the individual merging halos. 
As a substantial fraction of the Galactic GC population is old and metal-rich, it behooves a 
model of the GC population formation to describe the origins of these GCs. 



Discussion 



The AMRs constructed by iMarm-Franch et all (120091 ) and IPotter et al.l (I2OIOI . further 
updated in this paper) show that G alactic GCs split into two distinct branches at [Fe/H] 
> —1.5. iMarm-Franch et al.l ( 120091 ) interpreted this branching as representing two groups 
of GCs: an old group with uniformly old ages (a flat AMR) and a young group with a 
signiflcant trend toward younger ages at higher metallicities. They concluded that these 
groups likely originated from two diflFerent phases of Galaxy formation: an initial, relatively 



Marm-Franch et al. 



brief collapse and a second, prolonged episode driven by accretion. 
( 120091 ) found a small number of metal-poor ([M/H] < —1-5) G Cs with somewhat younger 
ages than the bulk of the metal-poor GCs. iDotter et al.l ( 120101 ) and the present study flnd 
all of the metal-poor GCs to be uniformly oldjfl but this diflFerence does not substantially 
influence the conclusion with respect to the proposed Galaxy formation scenario. 

A number of influential studies have used the HB morphology-metallicity diagram to 
infer the ages of stellar populations for which metallicit y information and Rood photonietry to 
at le ast the level of the HB a r e avai lable. For example, ICatelan fc de Freitas Pachecd ( 119931 ) 
and iLee. Demarque. fc ZinnI ( 119941 ) both used synthetic HB models to demonstrate that, 
under the assumption of constant mass loss, the Galactic GC HB morphology-metallicity 
diagram showed evidence for an age spread of several (^ 4) Gyr for a flxed He content and 
heavy element distribution. 



The use of the HB morphology-metallicity diagram must, however, remain a crude age 



^See Figure 9 of 



Potter et al, 



(|2Q1Q[ ) for a direct coraparison of the AMRs. 
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indicator at best for a variety of reasons. The lack of detailed abundance information, 
including the possible spread in He within a given GC, has important implications for HB 
morphology. The effect of a spread in the total He conte nt on HB morphology among 



different GCs has been explored by many authors (see, e.g., ICatelan fc de Freitas Pacheco 



D^Antona fc Caloi 



200^ 



I993I: iLee. Demaraue. fc Zinnlll994j ) as well as a spread in He within single GCs (see, e.g.. 



Another important factor in matching GC HB morphology to 
model predictions is statistical fluctuations. The ACS GC Treasury project photometric 
database contains GCs with anywhere from ^ 20 to several hundred HB sta rs, leading to 



uncer tainty in the median color of the HB of up to 30% in the sparsest GC (IDotter et al 
2OI0I . Table 1). 



Addressing the question of which para meters most s t rongly infl uence HB morphology 
using full CMDs for large samples of GCs, IPotter et al.l ( 120 lOl ) and ICratton et al.l ( 120 lOl ) 
reached the same conclusion: that age is the second parameter influencing HB morphology. 
The results presented in this paper lend further support to this conclusion. It is important 
to note that the term 'second parameter' used here applies specifically to the variation of 
HB morphology among different GCs , the context in which the term was originally used in 



the 1960's (e.g., Ivan den Berghl 119671 : ISandage fc Wildeyl 119671 ). 'Second parameter' is not 



intended to apply to the spread in HB stars within a single GC, for which th e same term 



i s ofte n em ployed, nor doe s this paper seek to address that important issue. iDotter et al 
( 120101 ) and ICratton et al.l ( 120101 ) both found evidence for at least a third parameter as 
well, suggesting that age and metallicity alone are not sufficient to fully characterize HB 
morphology. 



Figure [TT] shows the HB morphology-metallicity diagram of the full sample of 68 GCs. 



The le 
RGB ( 



't panel considers H B morphology as the median color difference between the HB and 



.^Potter et al.l 120101 ) a nd the right panel instead uses t he (B— R)/(B+V+R) metric of 

(ll989h as compiled by iMackey fc van den Berghl fcoosh . IPotter et all toid ) demon- 



strated that these two HB morphology metrics are strongly correlated, with the main dif- 
ference being that A(V — I) does not saturate as (B— R)/(B+V+R) does at ±1. In this 
context, it is clear tha t age is th e second parameter i nfluencing HB morphology as con- 
cluded bv lPotter etaD koid ) and ICratton et al.l koid ). 



While Figure [TT] is not well sampled over the full range of metallicity, it can still be 
seen that the transition from a red HB to a blue one happens over < 0.5 dex in [Fe/H] 
for a flxed age, or < 2 Gyr at flxed [Fe/H]. It is interesting to consider the ensemble in 
this manner, despite the fact that age uncertainties remain at or near the 10% level. It 
is a fortunate coincidence that we observe these GCs at a time when the HB morphology- 
metallicity diagram is rich with information. Turn the clock back ^ 2 Gyr and the HB 



Table 3. Results of isochrone fitting 



Cluster 


DMv 


E{B - V) 


[Fe/H] 


[tt/Fe] 


Age (Gyr) 


IC 4499 


17.07 


0.18 


-1.6 


0.2 


12.0 ±0.75 


Ruprecht 106 


17.12 


0.14 


-1.5 


0.0 


11.5 ±0.5 


NGC 6426 


17.74 


0.36 


-2.2 


0.4 


13.0 ± 1.5 


NGC 7006 


18.16 


0.04 


-1.5 


0.2 


12.25 ±0.75 


Pyxis 


18.64 


0.25 


-1.5 


0.2 


11.5 ± 1.0 


Palomar 15 


19.50 


0.38 


-2.0 


0.4 


13.0 ± 1.5 


Palomar 5 


16.86 


0.08 


-1.4 


0.2 


12.0 ± 1.0 
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Fig. 11. — A demonstration of the age effect in the HB morphology-metallicity diagram 
usi ng ages derived fro m the MSTO. The left panel shows the A(V — I) metric as compiled 
bv [Potter et al.l fcOlOh and this stu dy; the right panel shows (B— R)/(B+V+R) as compiled 
by iMackey fc van den Berghl ( 120051 ) . The shade of each point corresponds to its age with 
the oldest GCs shown as blue and the youngest shown as red. 
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morpholoKY-metallicit y diagram would have been populated only on the red side^as, for 
example, in the SMC ( IGlatt et al.ll2008l ). At such a time, the study by lSearle fc ZinnI ( 119781 ) 
would have revealed little variation in HB morphology among the halo GCs and, hence, 
little or no insight into the formatio n of the Galaxy. Only the truly peculiar GCs, such as 
NGC 2808 dPalessandro et al.lboiol ). miffht have shown a sizable variation. On the other 
hand, turn the clock forward ^ 2 Gyr and all but the most metal-rich Galactic GCs would 
have blue HBs. 



8. Summary 



HST/ACS photometry of 6 Galactic halo GCs, IC 4499, NGC 6426, NGC 7006, Palo- 
mar 15, Pyxis, and Ruprecht 106, were presented. The resulting CMDs were used to derive 
ages via isochrone fitting. The a^ e of Palomar 5 was derived usi uR the s a me se t of models 



and the combined photometry of ICrillmair fc SmithI ( 1200 ih and IStetson 



and metallicities of these 7 GCs were added to the AMR of Potter et al. 



(2000). The ages 
kom . The total 



number of homogeneously studied GCs in the sample is 68, excluding those 7 GCs imaged 
by the ACS Survey of Galactic GCs that are known to harbor multiple, distinct stellar 
populations. Divided at Rqc = 8 kpc, the inner Galaxy GCs exhibit a pattern of rapid 
chemical enrichment spanning two orders of magnitude in metallicity over a timescale of 1-2 
Gyr. The outer Galaxy GCs exhibit a much slower chemical enrichment that is evocative 
of dwarf galaxy cher aical evolution and i s reaso nably well matched by the semi-analytic GC 
formation models of iMuratov fc GnedinI (120101 ). 



Data presented in this paper were obtained from the Multimission Archive at the Space 
Telescope Science Institute (MAST). Support for this work (proposal GO- 11 586) was pro- 
vided by NASA through a grant from the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under NASA contract 
NAS5-26555. This research has made use of NASA's Astrophysics Data System Bibliographic 
Service as well as the SIMBAD database, operated at CDS, Strasbourg, France. 
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